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Abstract 

Analyzed is the effective potential of L>-dimensional thermodynamic Gross-Neveu 
model (defined by large-A leading order, 2 < D < A) under constant electromag- 
netic field (E ■ B = 0). The potential is derived from a thermal analogy of the 
worldline formalism. In the magnetic case, the potential is expressed in terms of 
a discretized momentum sum instead of a continuum momentum integral, and it 
reproduces known critical values of fi and T in the continuum limit. In the electric 
case, chiral symmetry is restored at finite T (with arbitrary fi) against instability of 
fermion vacuum. An electromagnetic duality is pointed out as well. Phase diagrams 
are obtained in both cases. 
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1 Introduction 



The Gross-Neveu (GN) model is known as the simplest model possessing the similar 
properties to QCD; dynamical breaking of (discrete) chiral symmetry, dimensional trans- 
mutation, asymptotic freedom and renormalizability {D < 4). Indeed, a number of GN 
models (and related models such as the NJL models) have often been studied as the proto- 
types of dynamical symmetry breaking in the various circumstances of finite temperature 
fl], 18, 19|, electromagnetic field and finite curvature [l^. The high tempera- 



ture behaviors of these models 0, |Tl|, [T^, ^} may potentially be interesting for exploring 
a dynamical symmetry breaking in (quark-gluon) plasma and in cosmological situations 
like early universe or a neutron star. 

In this paper, we examine the effective potential and chiral phase structures of the 
D-dimensional (2 < D < 4) GN model defined by the leading order of 1/A^-expansion in a 
magnetic B (or electric E) background field at finite temperature and density. In the case 
of D being slightly below 4, the plasma physics is an interesting application ||14| . In the case 



of D = 3 |n|, the physics of our model concerns planar condensed matter systems such as 
BCS theory: a magnetic field works as a catalyst for generating an energy gap in fermion 
spectra, through enhancing the interactions of fermions at small momenta even if the 
weakest (attractive) interaction. Gusynin et. al. suggested |I^ that this effect should 
be universal, i.e., model independent, in Z) = 3 and should be considered in the phenomena 
such as quantum Hall effects and high temperature superconductivity. Similarly in a 
supersymmetric model, constant magnetic field generates fermion's dynamical mass p^ . 

In the case of D = 2, magnetic field can not exist, however the GN model by itself cor- 
responds to a model of poly-acetylene. The first order phase transition of poly-acetylene 
from solitonic (broken) phase to metallic (symmetric) phase is formulated |^ as a func- 
tion of doping concentration (chemical potential) within the large- iV approximation. The 
critical value of chemical potential coincides with an experimental value, and recently a 
1/A^-correction to this result was obtained in ||^. Our study might explain some ther- 
modynamic behaviour of a poly- acetylene model under constant electric field in a finite 
temperature and density system. 
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Our crucial interest is therefore the question how universal the effect of magnetic 
catalyst on phase structures is in the setting of both finite temperature T and chemical 
potential /z in various dimensions. We study this question for an electric background 
field as well. In order to deal with these models all together, it is, of course, useful to 
analyze a dimensionally regulated effective potential. Such approaches were done in the 
cases of finite temperature |T^, curvature |T5[ and their mixture [jl^. We shall derive a 
universal expression of the effective potential as a function of D, T, fi and the external 



field parameter ^ (= ey l-B^ — where we have assumed B ■ E = 0) within the leading 
order of the 1/A^-expansion. We shall then present various interesting properties which 
can be derived from the universal effective potential. Although we have not completed 
our analyses yet, we confirm that our present results will be useful in order to establish a 
foundation for further developments in the various branches mentioned above. 

This paper is organized as follows. In Sect. 2, we present a simple and convenient 
derivation of the effective potential. There is no manifestation of exact fermion propagators 
of finite temperature field theory in the derivation. Namely, we do not use its Feynman 
rules to evaluate a fermion loop. Indeed, we apply an analogy of the worldline formalism 



p2[,|2^ to the finite temperature cases. Then we write down the proper-time integral 
representations for the critical equations which determine the phase structures in the ^- 
T-fi space. 

In Sect. 3, we discuss the cases of magnetic field dominant, i.e., ^s a real number. 
In the former half of this section, we analyze the zero-temperature limit of the effective 
potentials for finite ^ and fi. The potential is shown to be a discretized version of the one 
presented in [|^. Also, first order phase transitions are observed in large /z regions on the 
^-fi (T = 0) plane. In the latter half of the section, we discuss the phase structures. We 
find that a constant magnetic field works as the catalyst in 3 < -D < 4 at finite T and /z, 
and verify that the magnitude of catalyst effect becomes smaller as D or fi increasing. 

In Sect. 4, we present the electric dominant phase diagrams, i.e., the cases of a pure 
imaginary number. Analyzing the imaginary part of the effective potential, we show that 
chiral symmetry can be restored in high T situations (/i > 0) in 2 < D < 4, with vanishing 
imaginary part. We also show that some analytic results at T = /i = can be obtained 
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through the electromagnetic dual rotation ^ — » ip. 

Conclusions and discussions are in Sect. 5. In Appendices, some details on the ^ ^ 
limit are presented. One can find how to reproduce various known results at ^ = 



18, 20, 19i from our results in this continuum limit ^ — > 0. 



2 The effective potential in ^-T-fi space 

The model Lagrangian is 

£ = ^Pi^D^tP - - d^^, (2.1) 

where a is the auxiliary field and the number of fiavors. The fermions ip = {ipi}, {i = 
1, ■ ■ -yN), are coupled to the background gauge fields via the covariant derivative = 
dn — ieAfj,. The effective potential is given by 

V{a,O = V{a)-V{0) , (2.2) 

with 

Via) = ^or^ - In Betii^D^ - a). (2.3) 

This model does not possess a continuum chiral symmetry but a discrete one. After the 
effective potential is obtained as a function of the auxiliary field, the potential becomes 
the same form as that of the NJL model. 

We shall show below that the D-dimensional thermodynamic effective potential is given 

by 

where G2 is the elliptic theta function of second kind and tr[l] means the trace of gamma 
matrix unit (we do not need its explicit value). The D = 3 (magnetic dominant) case was 
derived in [|ll|]. Also, by letting T(= /?~^) and /i to be zero, this potential reproduces the 
potentials derived in the D = 3 magnetic dominant (real ^) cases and the electric 
dominant (pure imaginary ^) cases in D = 2 0, |25|, -D = 3 and D = 4 [Q. 
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2.1 Derivation of the effective potential 



In this subsection, we present a somewhat different derivation of the effective potential 
from that done in |]TT|, where they evaluated it in the Matsubara (imaginary time) method 
using a fermion propagator under a constant magnetic field. We shall derive the potential 
(|2.4|), modifying the worldline formalism |22|, E3[ into the thermodynamic case. 



Let us start with reviewing a few notes on V{a) for convenience of notation. In the 
Schwinger proper time method |2^, a basic expression for the potential (|2.3|) already exists 
in[| 

ds 1 „„2 



exp 



-ie.F,.a--iTrln(^^^ 
2 2 V esF 



(2.5) 



where tr and Tr mean the traces in the Gamma and Lorentz matrix spaces respectively. 
The trace of F^" is given by the following in respective dimensions: 2E^'^ in D = 2, 
2(^2 - ^2)" in D = 3 and 2[(p^)" + (p^ )'^] m D = 4, where 



pI 



±^{B^ - E^Y + A{E ■ B) 



If we assume E-B = in D = 4, and if the notation p = ey E"^ — B"^ is properly understood 
in each case (for example, B should be zero in D = 2), we can write those trace values in 
a unified form 



Tr 



[eF) 



2n 



fj,u 



2p 



2n 



(2.6) 



Similarly, the eigenvalue of ^ea ■ F can be expressed as ±ip multiplied by the degeneracy 
factor |tr[l] in each D, where we have again used E ■ B = in D = 4. In this way, we 
re- write 



tr exp 



1 



2 ^ 



Tr 



In 



sin(esF) 



esF 



tr[l]cosh(s^), 
sinh(s^) 



= 21n 



(2.7) 
(2.8) 



where i = ip = ey B'^ — E"^. Note that ^ can not be real in D = 2 because 5 = 0. Thereby 
the potential (U) becomes §, |, |, |, |n[ 

<-°° ds I 



s (47rs)^/2 



e-"' s^coth(sO. 



:2.9) 



Note that the external field dependence coth s^, the loop momentum contribution 
(47rs)~-^/2 and the mass term e are decoupled in this expression. 

Our proposal for introducing temperature and chemical potential is the replacement 
of the (normalization) factor {Atts)^^^"^ with a thermodynamic normalization. To see this, 
let us consider an alternative derivation of the second term on RHS of (|2.9| ). It is the same 
as the one-loop effective action of a spinor loop (with mass a) coupled to a background 
gauge (constant field strength) = ^x'^F^,^. In the worldline formalism, the effective 
action is organized in a form of closed path integral of one- dimensional Euclidean action 



22, 23 



T[A] = lnDet(i7^L>^ -a) 



X exp 



1 1-1 

dr^-x'^ + -^''ip + —iexFx — ietpFtp) 



(2.10) 



where x = x'^(r) and ip = iIj^{t) are the periodic and anti-periodic world-line (real) fields 
which are defined on a closed loop parametrized by r (0 < r < s), and x means drx{T) 
etc. The closed path integral amounts to the super determinant 



SDet-i/2 



^ + ^leFdr 







^^ 







^dr - leF 



Det(l - 2zeF(9-i)p 
Det(l - 2ieF(9-i)A 



"1/2 



where M is defined by the free part integral |22 



J\f = j) T>xT>ip exp 



1 2 1 -N 

vTa; +-il)i>)dT 

4 2 



and its explicit value is given by the formulae 



Pxexp 



x^dr 



4 



d^k 
(2^* 



T>ip exp 

Thus the effective action yields 



m = -^tr[l] / ^(47r.) 



-ip-^dr 
2 



= (47rs)-^/2^ 
-tr[l] . 

Det(l - 2ieFd-^)A 



Det(l - 2ieFd-^)p 



1/2 



(2.11) 



(2.12) 



(2.13) 
(2.14) 



(2.15) 
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The determinants can be evaluated on the Fourier bases {exp[27riA;r/s] |A; G Z, /c 7^ 0} for 
the periodic and {exp[27!'i{k + / s]\k G Z} for the anti-periodic functions |23|. Note that 
the zero mode in the periodic basis is excluded. We have 



2'7iik 



Det 



,sin(eFs) 



eFs 



(2.16) 



Det(l - 2zeF(9^^) = Det' fl (l - 2zeF 

fc=— 00,7^0 

and similarly 

00 

Detfl - 2ieFdZ^) = Det' TT f 1 - 2ieF — 

V ^ Ja ^li^V 27rz(fc + l/2) 

where Det' concerns the determinant on a Lorentz matrix. We then arrive at the expression 



Det'cos(eFs), (2.17) 



|23 



T[A] = -^tr[l] J ^(47rs)-^/'e-"'^Det'^[eFscot(eFs) 



(2.18) 



and this coincides with the second term of (|2.9| ) with the identification of (p.6|) . 

Now let us proceed to the thermodynamic case. We have to replace ko-integTal accord- 
ing to the Matsubara formalism 



TT 



ko -{2n + 1) - ifi, 



dkf) 1 ^ 



(2.19) 



Namely, we change the path-integral normalization for xq through discretizing fco-iiitegral 
in(|ll 



exp 



4"^° 



1 °° r 

- exp[-s{-(2n+ 1) 
1 2 2/i . 47r 



and (|2.13|) is then replaced by 



® Vx exp 









(4vr.)--^^e^-^e,(.|^,z.^; 



D-l 1 

T3' 



Since the zero temperature limit ^ 00 of RHS of ( p.2(]| ) is 

1 2 2/i . 47r 1 
hm -e^'^ 62 s-5-,2s— = ^= 



(2.20) 



(2.21) 



(2.22) 



( |2.13| ) is reproduced from ( |2.21| ). Equivalently, this can be verified by decoupling the zero 
temperature normalization part (47rs)~^/^, 

1 . 



Pxexp 



— / -X 

4 



ji Vx exp 1^— j -x^ 



+ 2(47rs)— ^(-l)"e 



/4s 



cosh(n/5/x) , (2.23) 



n=l 



where we have apphed the transformation 



e 



SfJ? 



oo 

= 1 + 2^(-l)"e-"'^'/^'cosh(n/3/i). (2.24) 



n=l 



Obviously, the second term in ( |2.23| ) vanishes in /? — >■ cx), and is ( |2.13| ) recovered. 



Let us remember that the thermal parameter jS ^ has nothing to do with the loop 
length s which is the Schwinger proper time of the spinor loop |2^. All integrations on 



internal momenta (fcj, /cq) are now done, and thus the operator determinants do not change 
any more, since dr acts on the proper time coordinate and F is a constant matrix. We 
therefore conclude 



1 /■rfs02(sf. 



IS 



4-7r ' 

m = --tr[l] J ^-j£;f^/'''-''''^^-otHsO. (2.25) 



Replacing the logarithmic part of (|2.3|) by this quantity, we get 

V..(-; f ) - + ^ r -^M^MOe-'^''-''-'- (2.26) 



2A 213 Jo s (47rs)(^-i)/2 
When D = 3, this thermodynamic potential coincides with the equation (B6) of Ref . |1TI| . 
Also, when ^ = 0, Vp^^{a] ^ = 0) coincides with the D-dimensional thermodynamic poten- 
tial discussed in |]TP[. When T = /i = 0, V^,^(c'";,^) is reduced to ( |2.9| ) because of ( ^.22| ). 
Finally, subtracting the zero point value of the potential, we obtain the announced result 

V{a- e, A /i) = Vp^,{a- - V^ft^(0; 0- (2.27) 

2.2 The critical equations 

In the following, we list the critical equations that are entire bases of our numeri- 
cal analyses of chiral phase structures. We should remember that the original effective 
potential ( |2.3| ) contains the ultra-violet divergence and also has an imaginary part for 
pure imaginary ^. When the imaginary part disappears, we shall impose the following 
renormalization condition (on the real part) 

1 



(2.28) 

(7=1 



and we can thus define the renormahzed couphng constant Xr in the following way 

ds 

= tr[lJKe / — 



(47rs)^/2 

The real part of renormahzed effective potential is then expressed in terms of as 



(2.29) 



ReVR{a;^,T,ix) 



1 9 1 r n 

+ -tr[l]Re 



2Xr 

1 , 



ds 



r<coth(sO 



ViTTs 2 / 2/i . 47r 



(47rs)^/' 

s— , is^l + a^e''{l - 2s) 



(2.30) 



The gap equation to determine the dynamical mass, i.e., the lowest energy configuration 
of 0", is then 

ds 

+ tr[lJKe / — 
\R ^ ^ Jo (4^ 

(2.31) 



1 ds 



X 



(47rs)^/2 

Q^{s^i,^s^^)+e-^il-2s) 



(3 ^ (3' (3^ 
where we have eliminated the trivial solution a = 0. 

Now let us derive the equations to determine phase diagrams in the C,-T-fi parameter 
space. The critical surface of first order is governed by the simultaneous equations 



ReVnia = m,^,/?,/i) = 0, 
d 



da 



Rel^R(a,e,/3,/i) 



0, 



(2.32) 
(2.33) 



where m denotes a non-trivial solution of ( p.31j ). Combining these equations, we obtain 



= Re/^ (^^^ecoth(.0 



ITS 



P 



2u 471 



) 



(2.34) 



One can decouple RHS of this equation into V{m; ) and a pure thermal part using ( p.24| ), 
however we stop writing further details, which will be discussed in sect. 3.1. The necessary 
condition to determine the critical surface of second order is generally expressed by 



lim AReV^^(a,e, /?,//) = 



(2.35) 



and this reads 




Re 



(47rs)^/2 



secoth«){-e^^ —^Q2{ 



+ e~ 



P 

2s) } + 2se~ 



2/i Att 



(2.36) 



Here, we have adopted the following value of the renormalized coupling constant 

which makes the chiral symmetry broken atT = yU = ,^ = with generating the dynamical 



mass cr = 1. Note that a solution of (|2.36|) does not always represent a true second order 
transition point if a first order transition occurs. We have to cut an irrelevant piece away 
from the true critical surface solving the equation of third order critical line. The third 
order critical line, which yields the boundary between the first and second order critical 
surfaces, is determined by the simultaneous equations ( |2.35| ) and 

2 



d_ 

'o^^'o \ da 



I^IttIJ ReVR(a,e,/3,/x) = 0. (2.38) 



The latter equation leads to 

One can perform any combination of limits and/or /i — > in the 2'nd and 3'rd 
order critical equations (|2.36|) and (|2.39|) . For instance, taking the limit ^ — in each 



equation, we obtain the following analytic equations respectively |T^, |19| (see Appendix 
A) 

/?^-r(l - ^) = -|(2vr)^-r(^^)ReC(3 - d\ + ^^), (2.40) 

and 

ReC(5-A^ + ^^) = 0, (2.41) 
Z Ztt 



where ( is the generalized zeta function. It is shown in [O] that these equations reproduce 



the critical equations derived in [0]. The limit /5 — > cxd can also be taken by applying 
( |2.22| ) to ( |2.36| ) and ( |2.39| ), and their limits are valid unless any I'st order phase transition 
takes place at T = 0. In contrast with these cases, the effective potential on the ^-/x 
plane (and hence the I'st order critical equation there) can not be obtained for fj^ ^ 
by naively applying ( p.22| ). In order to reproduce a precise limit, we have to collect all 



contributions from under the Fermi surface at T = 0. This can be achieved by an insertion 
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of the step function 6{jj, ~ k) in a loop momentum integration at a very beginning of 
computation. This prescription is, however, unclear in our present formulation, since we 
already integrated the loop momenta. We observe a method how to deal with this problem 
in the next section. 



3 The magnetic dominant cases 

In this section, we discuss chiral phase structures in the magnetic dominant cases for 
finite T and fi. We formally include D = 2 case also, because it will be useful when ^ is 
rotated to ip and because it is convenient to see a D-dependence of magnetic catalyst in 
the region 2 < D < A. Since C, is real, the effective potential does not possess an imaginary 
part, i.e.. 

Re VR{a; e, /?, /x) = Vn{a; e, /5, /i), (3.1) 

and there is no more intricacy of critical equation analysis (up to the T —>■ limit). In 
T^, the D = 3 case (/i = 0, T 7^ 0) with a constant magnetic field was analyzed, and the 



magnetic field turns out to play a role of strong catalyst of chiral symmetry breaking. We 
shall observe this feature in the ^-T-/i space for various values of D. Each critical surface 
is determined by numerically analyzing the I'st, 2'nd and 3'rd order critical equations 
presented in the previous section. As to the zero temperature limit, expected from [Q, we 
shall show that first order phase transitions exist on the ^-fi plane at large values of fi. 



3.1 The effective potential on ^-fi plane 

First we study the zero temperature limit for finite ^ and fi. As mentioned at the end 
of sect. 2, the special care is necessary to analyze the effective potential (not gap equation) 
on the ^-fi plane. In addition, we have to remember that C,~^ is relevant to the magnetic 
length eB (or C, the discrete momentum unit), which is known as the radius of cyclotron 
motion of an electron. This assures the magnetic lattice structure, although it can not be 
seen in the theta function representations of the effective potential (|2.3CI| ) etc. 

The effective potential ( p.27| ) with the renormalization ( p.29| ) can be written as the 
sum of the zero-temperature potential (|2.2| ) and the pure thermodynamic remainder V by 
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means of the formula ( ^.241 ): 



V{a■i,P,^^) = V{a■0 + V{a■i,P,^^) 



or (q.v. (|23) and (^^) 



Vp^,{<j-i) = V{a) + Vp^^{a-i), 



where 



V{a) 



tr[l] ds 



2Xr 

V>,^(a;0 = tr[l] 



2 Jo (47rs)^/2 
°° ds 



coth(s^) 



a'e~%l-2s) 



/o (47rs)^/2 
Using the following formulae 



ecoth(so-e-^'^ E(-ire 



cosh(n/3/i) 



n=l 



oo 



coth(sO = 1 + 2^6 
1=1 



-2ls( 



dss" exp 



7S 



2 g) ' /C(2v/^ 



we can perform the s-integrations in ( p.4|) and (|3.5| ) just similarly to the D = 3 
The results are 



a- . tr[l]e 



2A 



+ 



R 



2(20--^C(l - y, 1 + ^)r(l - + a^-^l 



2(47r)^/2 

+2a^(20--^C(2 - f , 1 + ^)r(2 - ^) + a^r(2 - ^) 
-4a2(20^-^C(3 - f , 1 + ^)r(3 - - 2a^r(3 - 



and 



tr[l]e 



^ 1=1 



(47r)^/2 

where we have introduced the compact notation 

C»;3(a) = 2/ rfss-^/2E(-l)"cosh(n/3^)exp 



n=l 



^ 4s ' 



4E(-l)"cosh(n/?/i) 



?i=i 



2a 



l-D/2 



KD/2-i{npa) . 



11 



Although we have performed the integrations w.r.t. the proper time s, it is rather conve- 
nient to consider an integral representation in order to observe the low temperature limit 
of 13(a). Using the formula 



^^^^^"r(z. + i) A 



e-'\r -ly-^dt , Rez/> — , Rez>0 



(3.11) 



we get the following integral representation 



-cr 



D-2 



CsO poo 

^(-)"cosh(n/?/i) / e-"^'^*(t'-l)^(it . (3.12) 
— 1 1 



V 2 n=l 

Performing the summation over n and changing the integration variable t = 1 + 1' /a, we 
obtain 



2^/? r r (^^ + 2c^i) 



r(^) 



-D-3 
2 







1 + e'^(*+'^+'") 



(it + (/i — ^ yU) 



Ofj{0) is given by (a derivation is in Appendix A) 



2 /27r\^-^ /3-g 
^Ui I 2 



)ReC(3-D,^ + 



t— 



(3.13) 



(3.14) 



Here, we stress that (|3.9| ) is a discrete version of the ^ = potential obtained in |T9| 



\/(a;0,A/i) 



tr[l]2y^ 



(47r)^/2r(^)/3 Jo 



:i + e- 



-/3(Ai+Va2+fc2) 



)(l + e- 



'I + e-'3(''+'=))(l + e-z^e^-/^)) 



(3.15) 



This can be easily understood in the case of D = 3, since Op{a) takes the following simple 
form 



0^(a) = -2^1n 



(l + e-/3(a+M))(i+g-/3(<x-/.)^ 



(3.16) 



and hence we have (see also [11, E 



tr[l]^2v/?rl (l + e'3('^-'"))(l + e~'^('^+'^)) 
' (47r)3/2/5 l2 (l + e/3/^)(l + e-/3/^) 



(1 + g/3(A«-V<^2+2i^)^^^ _^ g-/3(A«+V'^'+2«C)) 



(3.17) 



This is exactly a discrete summation form of ( p.l5| ), and the correspondence to ( |3.15|) is 
understood by the following translation rule from discrete spectrum to continuum one 



2/^ 



12 



oo 



(3.18) 



kdk . 



Remember, in the continuum expression, we have to (i) insert a step function 6{fi — k) 
into the above momentum integration and have to (ii) keep fi < a in order to observe a 
correct potential behavior in (3 —>■ oo with finite fi; namely 

tr[l]2v/^ 



lim V{a; 0, /3, /i) = ^n/2r(D-i^ I 



dkk^'-'iij - k) . 



(3.19) 



Due to the correspondence (|3.18|) , we may similarly assume an upper bound on the discrete 
momentum sum; I < L = [^i^ /2^]g, where [ ]g is the Gauss symbol. 



lim V{a] ^, (3, fx) 

p— >oo 



tr[l]e20F 



(47r)3/2 



i/i + 5:(/i-V^) 



(3.20) 



This potential is still clearly a discrete form of ( |3.19D for D = 3 under the identification 

of ( CT) . 



Let us return to the generic D case. The above cutting procedure is also valid for the 
generic D case, because we naturally have a step function in 



Ooo(a) = lim 0(3{a) 
20F 



(3.21) 



-^(2a)^(/i-a)^F(^-,^-;^-;^)0(/i-a) , 



where we have applied the integration formula ( [B.3D . Owing to the presence of the step 
function 6'(/i — a), the upper bound L is naturally introduced here, and we figure out 



lim V[a] (3, fi) 

p— >oo 



'"^^^^ riOoo(0)+EOoo(V^) 

^ 1=1 



(3.22) 



(47r)^/2 

This is the discrete form corresponding to ( |3.19|) . The explicit forms of (9oo(cr) in each D 



are 



-p3Y-0(/i-a) X < 



arccosh (/i/ci) 
/i — cr 

l/i-y/yU^ — cr'-^ — |cr^arccosh (yu/ a) 



for = 2 

for D = 3 (3.23) 

for D = A, 
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and 



Ooo(O) 



2^ jj 



D-2 



D 



(3.24) 



The Ooo(O) is singular when D = 2, however it vanishes under the 'physical' limit ^ — > 0. 
In Appendix B, we show how the discrete form ( ^.22| ) reproduces the continuum one ( pj.l9D 
in the limit ^ -> via ([3181 ). 

Now let us examine the phase structure on the ^-/i plane. It is enough to analyze the 
I'st order critical equation defined by (2.32) and ( [2.33 ). The reason will become clear 
at the bottom of this paragraph. Since the pure thermal part ( |3.22|) does not depend on 
a any more, the gap equation is the same as the one for the zero temperature potential 



2s) 



(3.25) 



(Of course, this can be derived from ( p.31|) with the use of (|2.22| ).) Eliminating the 
renormalized coupling from ( |3.4| ) with using ( p.25| ) , and applying the following formula Q 



dss'^-^e-'^' coth s = r(/i) f2i-^C(/^, 1 + -) + P'"] , Re > 1 
1-2-1 



(3.26) 



V{(t; ^)\(j=rn turns out to be 
y(m;0 - 



tr[l]2^/^e^-^ D 
2(47r)^/2 ^ 2 ' 



-ir?{\ C(2 ,1 + — 

2 ^ 2^ V 2 2^ 



+ ec(i--,i + 



D 
~2 



Combining this and (|3.22| ), we obtain the I'st order critical equation on the ^-/x plane 

-2 



(3.27) 







rri 



D 



1=1 



+ m2(2-D)c(2-Z},l + ^) + ^y^(20'"^m 
, 4«(l4.1.f)-4«(l4 



D-2 



(3.28) 



where m should be determined from ( p.25| ). (For example, we can take m = 1 as a solution 
if we choose, instead of (|2.37| ), 



tr[l] 



°° ds 



(47rs)^/' 



coth(s^) ■ 2se 



(3.29) 



This formula is derived from the expansion (|3.6D and the definition of zeta function 
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This coupling choice makes our critical equations much simpler, however it is not essential 
for the numerical analyses of phase structures discussed later). Under the choice ( p.37| ), 
the mass Q is determined by 

= 2-e(3-D)-em^-^-(20^-^C(2-f,l + ^) 

+ (20^-^C(2 - Y> 1 + ^) - (20^-^(4 - D)C{3 - f , 1 + (3.30) 



The dynamical mass determined by ( p.30| ) grows as the strength ^ of the external magnetic 
field becoming large (see Fig.l), whose feature coincides with the one argued in 0]. Note 
that we can not take the limit m —>■ without circumventing divergence for D < A in this 
equation. This means there is no 2'nd order phase transition on the ^-fi plane. 

Substituting the solutions of ( |3.30| ) into ( |3.28| ), we find I'st order transitions (see Fig.2) 
as the solution of (|3.28| ). The equation ( p.28|) can be shown to coincide with the following 
equation |]19| in the limit ^ (The proof is in Appendix C), 

which has the solutions /i = m/\/2 for D = 2 and jj, = m for D = ?> [113, BO]. 



(3.31) 



m 




Fig. 1: The dynamical mass m for 
magnetic dominant cases on T = 
plane. 




0.0 



5.0 



10.0 



Fig. 2: The values of critical chemical 
potential yUc on the T = plane clS 3j 
function of the external magnetic field 



^For example, a solution is m = 1.02107 for ^ = 0.3 in D = 2>. 
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3.2 The phase structures 

In Fig. 3, we illustrate several critical curves which represent various sections of the 
second order critical surface in the C,-T-fi space for D = 3. In this case, the point T, /i) = 
(0,0,1) is the unique tricritical point, since (|2.41|) has the solution {j3fi)~^ = 0, which 
means T = (for /i 7^ 0). Fig.4 is the case of D = A — e, where we take e = 0.5. No 
tricritical point appears in this case. We could not let e to be much smaller than this 
value due to a singular situation of numerical computations. However we can easily read 
some signs of the magnetic catalyst and the disappearance of the tricritical point, if we 
extrapolate from the following behaviors of larger values e = 1.5 and 2.0. 

First, let us observe the D = 2.5 (e = 1.5) diagrams in Fig. 5, where the tricritical 
point starts appearing around 0.75 < fi < 0.76. These tricritical points are determined 
from intersections between two curves solved from ( p.36[ ) and ( |2.39| ) like the way shown 
in Fig. 6. Next, let us watch the D = 2 diagrams (Fig. 7). The tricritical line (points) and 
first order transitions are observed for 0.61 < fi. The phase diagram is very similar to 
Fig. 5, and these diagrams exhibit that strong ^ makes phase broken (magnetic catalyst) 
even in these unphysical dimensions 2 < D < 3. Now we gather the tricritical points of 
each dimension in Fig. 8. It clearly shows that the tricritical line moves toward a lower 
temperature region as D increasing. And D = 3 is the terminal case that the tricritical 
point appears. We therefore deduce that there is no tricritical point in D = A — e for small 
e whilst observing the magnetic catalyst in finite T and fi. 

Let us summarize the phase structures. We have observed that the critical surfaces in 
D = 3 and 4 — e are of second order, and magnetic fields tend to make the chiral symmetry 
broken even at finite T and fi. The large values of fi tend to revive the chiral symmetry, 
and the critical /i grows as ^ increasing (for all T). There are first order transitions not 
at finite T but at zero temperature shown in Fig. 2, where the ^ = axis is singular and 
discontinuous. 
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Fig. 3. The critical lines for magnetic 
dominant cases in D = 3. 




Fig. 4: The critical lines for magnetic 
dominant cases in D = 4 — e, where e 
is chosen to be 0.5. 




Broken 



1.0 2.0 3.0 



Fig. 5: The critical lines for magnetic 
dominant cases in D = 2.5. 
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T 



0.5 



0.0 



Symmetric 



C3 



/ C1 Broken 



0.0 



1.0 



2.0 



Fig. 6: How to determine the tricriti- 
cal point A. This example is the case of 
fi = 0.7, D = 2 (magnetic dominant). 
The C2 and C3 are the solutions of eqs. 
(|2.36|) and (|2.39|) . The CI represents 
the solution of ( p.34|) . 
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Fig. 7: The critical lines in D = 2 for 
magnetic-like external field. The T = 
plane is singular. 
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Fig. 8: The tricritical points projected 
on the T — ^ plane for magnetic domi- 
nant cases. 



4 The electric dominant cases 

4.1 The imaginary part of effective potential 

In contrast with the S-dominant cases, the effective potential ( p.4|) gives rise to an 
imaginary part due to the instability of fermion vacuum in the i?-dominant cases ^ . 

In fact, the function s^coth(s^) in the integrand of (|2.4|) is now spcot(sp) with changing 
^ = ip, and the integrand then contains an infinite number of the singular points s = 
nir/p, {n = 1,2, ■ ■ ■) accordingly. Then we must displace the contour of integration (— oo < 
s < oo) by ie along the imaginary axis of the complex s-plane. The contribution of each 



singular point is taken into account by making use of the ie rule (for example, |25]) 



1 P J mi\ , , 

mb\s , (4.1) 



p p \ K 



where P stands for the principal-value symbol. By virtue of this rule, the imaginary part 
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of the effective potential consists of 

p{D-l)/2 oo 



n=l 



oo 



-tr[l] 



^^Re V exp {a^ 



pf3 p(3^ 
2m + 1 



n- 



X 



»( 



D 



-1 2 



m=0 

, , , ,2m + 1 



P 



(4.2) 



where ^{z;s,a) is the Lerch transcendental function defined by ( p.7|) . This can be re- 
garded as a probability, per unit space-time volume, for a fermion pair-production from 
the vacuum by an external E field. Taking the T = /i = limit 



Mm lmVp^f,{a; p) 

1 ,lJ,—fO 



V2' 2p' 



2^+%^-i 
we reproduce the known results p|, 

tr[l]p/(87r) ■ ln(l - e^'^'^^P) for = 2 

-tr[l]p3/2/(47r)2 . e-^'^'/^<l>(|; ig, 1) for D = 3. 
in the following form using the Lerch transformation ( p.9|) . 



(4.3) 



lmV{a 

One can rewrite 

Iml^(cr 



(4.4) 



4(27r)^/2 



-|- exp 



2 ' 

-.(1 - ^) 

2 ^ 2^ 



rvr 



exp 



^(1 



2 ' 



C(l ,i— ) 

2 ' 2p^ 



TT 



CO -fa 



^ — 



(4.5) 



We shall search for the lowest energy configuration of a solving the critical surface equa- 
tions derived in sect. 2. This ground state solution cr(i?) minimizes the real part of the 
effective potential, and the chiral symmetry is broken accordingly. On the contrary, we 
expect that a strong external field might induce a restoration of chiral symmetry. However, 
we can not conclude immediately whether or not the chiral symmetry is restored even if 
the fermion mass vanishes, because the effective potential possesses the imaginary part 
that represents a decay of "symmetry vacuum" to a fermion pair. Indeed, in the case of 
T = /i = 0, the chiral mass term a{E){ipilj) does not vanish even when (j{E) = in D = 2 

i- 

Hence we determine the symmetry restoration by the criterion whether the chiral mass 
term vanishes or not 0, ^ 

lima(p)Im\/'(a(p);p,/?,/i) = 0, (4.6) 

p— >Pc 
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where we have apphed the fact that {ipip) is related to the Im V by 

{ip^)a=a{p) = iNlmV'{a{p);p,p,fi) - yo-(p), 



(4.7) 



and (t(p) denotes a solution of gap equation, i.e., ReV = 0. To check this criterion, it 
is sufficient to see whether or not the Im V'{a; p,P,p) is finite as cr ^ 0. From (^]2|), we 
have 



lmV\a;p,p,fx) = ati[l] 



The upper bound for the summation value in ( |4.8D can be estimated in the following way 



3-D _rnv/2 ,,2 



62 ( 



2mTp Amr' 
i- 



p(3 ' p(3^ 



71=1 

E°° 3^ ^ r TT 2 / 2m + 1 . 2t 
n 2 ^ exp^-n-|(T +(7r — ip) } 



71=1 
("OO 



7?1= — 00 

00 



, 3^ ^ r 2 , 2m + l . ,2, 
< / dxx 2 exp — X— |cr + (tt ip) | 



7Tl= — 00 



Performing the x-integration 



RHS (4.9) = Ti 



we find it finite as cr ^ 



5-D, ( p 



vr 



i=— oo 



ItoRHS (4.9) = 2r(^) (^^^ 



g-5 

2\ 



The criterion is therefore satisfied 



(4.9) 



(4.10) 



(4.11) 



lim aim V'{a] p, /?, p) = 0, 



(4.12) 



and we conclude that the chiral symmetry is restored for finite T and p in 2 < D < 4. 

Note that RHS of ( 4.11 ) vanishes right at a tricritical point, whose necessary condition 
is (|2.41|) . This is not the case with 2'nd order transitions and suggests that the multi- 
criticality can be mightier than the vacuum instability. As a by-product, we obtain a 
necessary condition for the tri-criticality on the imaginary part 



lim — Im V'{a] p, (3, p) = 0, 



(4.13) 
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and it leads to 



1 3^ /2nna 4n7r^ 



0. 



(4.14) 



It is interesting to note that the powers of a in front of Im V in ( [4. 12|) and (|4.13|) differ 
by two from each other. Hence ( |4.12|) and ( fl.l3| ) may be called a kind of discrete versions 

of ( p:^ and (p:^ ). 

In the zero temperature cases, one may derive exact expressions for Iml^' from ( [4.2|) 
by taking the T,fi — > limit like done in ( [4.4| ). Here, we only put a brief statement about 
their power dependence on a in this limit. Using 



1 ^ r ^ r 2 / 2m+ 1,2V 
i™ ^ exp -x-{a + (vr — - — ) } 



e p 



V4 



vra; 



we obtain 



hm IrHS (4.9) = ^ ^ 



2 ' 



and its power dependence on a exactly coincides with known results 



(4.15) 



(4.16) 



lim IvdV' {a; p, (3, ii) ~ < 
r,/i->o 



const. 



for D = 2, 
for D = 3, 
for D = 4 - e. 



(4.17) 



Namely, chiral symmetry is restored in D = 3 and 4 — e, and broken in D = 2 



4.2 The duality and phase structures 

In this subsection, we present phase diagrams with analyzing the real parts of the 
critical equations ( p.36| ) and ( p.39| ). Also, we point out a dual relation between the electric 



and magnetic dominant effective potentials (at T = 0) under the rotation ^ —>■ ip. We 
follow the same method employed in |] in order to evaluate the real parts (principal 
values) . 

Let us introduce the following compact notations for the real part effective potential 
(12.301 ) (to avoid lengthy expressions). 



ReVn{a;p,P,p) = ^a' + ^^^^{hia)-hiO)) , (4.18) 
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where 



and 



/i((t) = P/ -^spcotsp- f{s,(3,fi) 



s 



D/2 



(4.19) 



.(.-so 



2 „,,2v/47rs 



62 (. 



2/i . 47r 



(4.20) 



Furthermore, spht /i(cr) = A{a) + -B(cr) into the following two pieces for convenience: 



A{a) 
B{a) 



D/2-1 



dr 



(47r)«/2 



r 



l-D/2 



r 



cotr ■ f{-,p,p) 



(4.21) 



oo /•n7r+7r/2 
D/2-1 ^ p / 



_D/2-l 



„tl ^n.-V2 (47r)^/2 



2ri-^/2cotr-/(-,/?,/x), 



E 



■ cot r 



r 



+ n7r)^/2 



- - (r -1 



(4.22) 



„^,7o (4vr)«/2 

These quantities A and B have another expressions. Namely, the integration in ( [4.21] ) can 
be performed 



A{a) 



2k 



-2EC(2fc)(^ 

1- D 7r,2 .2m+l 



P m=0 



/5 



X7f- h 2A;, — |cr^ + (vr- ^ 

+ a^(l-- + 2^,-)-2.%(2-- + 2^,-) 



(4.23) 



Here we have applied the formula 

r'r/2 

dxx^e~°'^ cot X 



fc=0 



TT 



(4.24) 



where C,{z) is the zeta function, and 'y{z,p) the incomplete gamma function defined by 

rP 



'y{z,p) = I e~H'-'dt Rez>0 . 
Jo 

The summation of the second line in ( ^.22| ) can be performed 



(4.25) 



n/2 



dr cot T ■ 



2Re ^ exp 



m=0 



r + TT , 2 . 2m + 1 . xo T 



P 



TT 



r ., , 2m +1 . X 9 . T 



TT 



TT 



-{-+-)/p$(:^^_2;— ,l + -)-(r 
^2 '2p' ^ 



-r 



TT 



(4.26) 
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Similarly, we rewrite the second and third order critical equations ( |2.36| ) and ( p^.39D in a 
straightforward way, 



/2(a) + 2r(2-:^) = 0, 



IM = 0, 



(4.27) 
(4.28) 



where the explicit forms of I2 and I3 are given in Appendix D. 

For simplicity, let us begin with analyzing the T = fi = case. In this case, ( [4.20| ) 
simplifies because of (|2.22|) (If one wants to consider the /U 7^ situation, it is necessary 
to follow the similar procedure analyzed in sect. 3, and see also [^). Let us write down 
Ji and l2- 



r-K/2 

( ' drcotr- 



2 zp 71 

+aVe-(^+-)/^$(^-l;f,l + -) 

Z zp 71 



P 

D T. 



a 



P 

p 2 / dr cot T ■ 
Jo 



2p' 



TT 



-o -D ( r , r , V 

r- + >r-{C(y.l + -)-C(^,l--)} 



and differentiating / 



d 

P 



cr=0 



/ drcotr- 



^i-^/2|_i + e-^/''(l--)| 

_,,i-^/2|_^(^ -1,1 + -)+ e-(^+-)/''$(^ - 1; f , 1 + - 

^ ^ 2 7\' ^ Z zp 71 



p V2 '2p' TT 



■)} 



We obtain the following critical p-values as the solutions of (|4.27|) with ( [4.31|) 



(4.29) 



(4.30) 



(4.31) 



0.583 


for 


D 


= 2 




< 1.238 


for 


D 


= 3 


(4.32) 


3.394 


for 


D 


= 3.5 , 
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and we verified that the real part of the potential given by ( |4.18| ) with ( [4. 291 ) and ( [4. 301 ) 
behaves consistently. Here we add a few remarks on these results. From D = 2 to 
D = 3.3, whose pc = 1.965, all p^-values excellently coincide with the solutions of the 
following analytic equation 







2-(2p)T-ic(2--)cos[-(D-2) 



(4.33) 



+Re 



(2v)*-'C(2 - f ■ 1 - ^) + (2>P)*-^(f - 4)C(3 - f , 1 - ^) 



2' 2p^ ' ' 2 

and this still exhibits a good agreement | at D = 3.5 {pc = 3.395). We found this equation 
from ( |3.30| ) through rotating ^ ip, taking real part and the limit m 0. The divergence 
stemming from m now drops off since it becomes a pure imaginary number by the 
rotation. Similarly, we verify that those critical values are consistent with the following 
potential derived by the same rotation procedure from (|3.8|); 



V{a;p) 



a 



2A 



+ Re 



R 



tr[l]ip 



2(47r 



)Z)/2 



D 



2a^(2^p)T-\(2--,l--)r(2-- 



2p^ 



D 



4a\2zp)--X{3 - -, 1 - -)r(3 - - 



2p' 



+ 2(2^p)f-ir(i-^){c(i 



2 



(4.34) 



2 2p' 2 

Note that this rotation procedure is not always allowed in general cases (for example in 
the thermodynamic part V)^ because we have used the expansion ( |3.ti| ) (or the formula 
(|3.26| ) derived from (|3.6| )), which is not valid for a complex analytic function. In this sense, 
the above coincidences are non-trivial results (although they are rather trivial before the 
expansion is applied). 

Now we discuss the phase structures in the p-T-p space. In Fig. 9, we depict phase 
diagrams for several p values in the case of D = 2. As mentioned above, the vacuum state 
is not invariant at T = p = 0. In this sense, the critical value pc divides a phase between 
massive states (broken phase) and quasi- st at ional states (quasi-symmetric phase), and the 
total probabihty of decay per volume is given by B 



r = — lim 2alm V'((j; p, oo, 0) = — . 



(4.35) 



^The deviation of pc grows as D increasing in the region D > 3.5. We believe this is related to a 
singular problem in numerical computations like encountered in sect. 3.2. 
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The similar situation is expected for finite /i (with T = 0) as well, and the F changes into 
r + AF, where AF is formally given by 



AF = - lim 2alm V'{a; p, oo, fx) , (4.36) 



although we have not estimated its explicit value yet. On the other hand, as discussed 
previously, the chiral symmetry can be restored (i.e., F = 0) when both of T and /i are 
finite. 

In Figs. 10 and 11, the cases of D = 3 and 4 — e (e = 0.5) are illustrated. As shown 
in ( [4.17|) , these cases have a stable (chiral symmetric) vacuum at T = /i = 0. Also, in 
the situation of finite T and /i, the chiral symmetry is restored similarly to the D = 2 
case. In each dimension, the critical surface goes down to a low temperature region as fi 
increasing. We could not find any tricritical point within a region of significant numerical 
data. In particular the computation did not work well below around T = 0.1, and we 
hence extrapolate the critical lines in Figs.9-11. 
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Fig. 11: The critical lines for electric 
dominant cases in D = 4 — e, where e 
is chosen to be 0.5. 



5 Conclusion and discussion 

In this paper, we found various interesting properties of the thermodynamic effective 
potential of the large-A^ Gross-Neveu model coupled to constant electromagnetic fields 
perpendicular to each other: worldline approach, magnetic discretization, electromagnetic 
duality, critical equations and chiral phase structures. 

First, we summarize the analytic properties of the effective potentials. In sect. 2, we 
presented a new derivation of the effective potential from the point of the worldline ap- 
proach developed in recent years W\ The similar formulation had been applied to 



constant curvature cases in |Tj], but we could not clarify the correspondence to other 
calculations. However, we emphasize that it is interesting that our method proposed in 
sect. 2 leads to the exact results |TT], Another spectacular of our potential can be seen 
in sect. 3. Our potential reveals the role of ^ as a discrete momentum unit even in the 
situations under finite temperature with keeping /i finite, and in particular, the discrete 
potential of D = 3 completely resembles to the continuum potential In other words, 
our D-dimensional thermodynamic effective potential is composed of the sum over discrete 
energy spectra e„ = zti/m^ + 2ri^, n E Z. We also argued how to take the /5 — > cxd limit 



26 



of the potential with ^ and /x kept finite, and we described the precise continuum hmit 
^ — s> 0, which gives rise to the known resuhs of fic- It means /i is properly taken into 
account in our derivation. 

We found in sect. 4 the tri-criticality condition for the imaginary part of the electric 
dominant effective potential, and also found the analytic expressions ( [4.33|) and ( [4.34|) 
which are dual to the magnetic dominant effective potential and its gap equation at T = 
fi = 0. When the effective potential possesses the imaginary part, the vanishing conditions 
on the first and second derivatives of ReV^ are not sufficient for defining second order critical 
and tricritical points. It is necessary for second order critical points to satisfy the condition 
( 4.12| ), which is a vanishing condition on the chiral mass term crlmV^'. Similarly, ( 4.13| ), 
cr~^ImV^' = 0, should be satisfied at a tricritical point as well. These vanishing conditions 
on ImV may be regarded as a discrete version of the real part conditions ( p.35| ) and ( p.38| ). 

Based on all these properties, we investigated the phase structures in various cases. 
First, in sect. 3, we showed the phase diagrams for magnetic dominant cases, and observed 
the strong catalyst effects of magnetic fields. When P = oo{fi^0mD = 3 and 3.5), the 
symmetry restorations are not the second but the first order. In these cases, the chiral 
symmetry becomes broken in a strong magnetic field region. This property is the same 
for finite f3 (with /i 7^ 0) also. These features do not change if we formally extend to the 
2 < D < 3 regions, from which we actually got some insights toward a region D being 
near to 4. As pointed out in ^8\, the effective potential of D = 2 by itself resembles to an 
electrostatic potential problem of the charge distribution s coth s along a one-dimensional 
line. It means that such an electric distribution rather acts as a "magnetic" catalyst, if 
we could realize the distribution in a one-dimensional laboratory. 

The phase diagram analyses for the electric dominant cases are more complicated. 
First we have to test whether the imaginary part of the effective potential vanishes as a 



goes to zero (in other words, the chiral mass term cr{^jjip) vanishes or not, q.v. (|4.6| )). At 
T = /i = 0, this can be easily checked by integrating ( |4.17] ) w.r.t. a, and it shows that 
the lmV{a) behaves like a^"^. For 2 < D < 4, the chiral symmetry gets restored. For 
D = 2, the ImV^' diverges as cr — 0, and the imaginary part ImV" never vanishes 0. 
However, at finite (3 and fi, we showed that the chiral symmetry is restored in 2 < 
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D < A based on the same criterion. This suggests that a symmetrization effect by high 
temperature can be stronger than the vacuum instabihty by an external field. As to p-fi 
plane, we have not yet succeeded in formulating the T = limit with /i being kept finite 
(In this sense, the achievement of the magnetic case was remarkable). When the fermion 
density is large, a density effect might enhance the chiral symmetry, and consequently, 
the instability would be suppressed. In order to address this question (at T = 0), it 
is necessary to study the I'st order critical equation like done in ( 3.28| ). An interesting 
problem is whether or not a dual relation like the one between ( p.30|) and (|4.33| ) holds in 
this case, because it contains a non-trivial part, the thermodynamic quantity V, which 
has non-analytic dependence on /x. 



A A proof of (3.14) 



First, we show how ( |2.36| ) and ( |2.39| ) reproduce ( |2.4CI| ) and ( p.41| ) respectively. Then 
we present a proof of (|3.14|) . In the limit 0, (|2.36|) becomes ( |2.40|) in the following 
way 
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and similarly ( |2.39| ) reads 
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Let us turn to the proof of ( p.l4| ) making use of the above calculation. Recalling the 



28 



definition of by ( |3.1(J| ), and using ( ^.241 ), we liave 

Tlie top line of ( |A.1| ) can also be written in the form 
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li 2 < D < 4, the following integral is convergent 
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and therefore (|A.3|) amounts to 
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Comparing the RHS of this equation with the bottom line of ( |A.1| ), we have proven ( p.l4| ) 



B The ^ ^ limit of (3.22) 



We prove that (|3.22|) is reduced to ( p.l9| ) in the limit ^ ^ 0. Obviously, the I'st term 
(9oo(0) in ( p.22| ) drops in ^ 0, since Ooo(O) possesses no dependence on C,. Hence, we 
have only to prove the following equality 
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Under the correspondence (3.18|), LHS of the above equation can be transformed into 
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and applying the formula 3.194.1 |29 
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we get 
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LHS (5.1) = I I't^-^f. - trF{^, 2; 3; -2^ 



29 



Now changing the integration variable, t = 2fi/{w + 2), 

LHS (5.1) =/i^2^-3 r dww''{w + 2)-'-''F(^—^ 
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and applying the formula 7.512.10 ||29 
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(where Re7 >0, Rea — 7 + (T>0, Re/? — 7 + cr>0, |argz| < vr), we find 
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Using a Kummer transformation to make the argument of hypergeometric function from 
— 1 to 1/2 and the following formula 
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we have 
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Finally substituting this into ( [B.7| ), the proof ends 

ir(3)r(D)r( 



LHS (5.1) =/x^^ 



fx 



D 



4r(D + i)r(^) D(D-i) 



(B.8) 



:b.9) 



(B.lO) 



C The ^ ^ limit of (3.28) 



We prove ( |3.28| ) goes to ( p. 311 ) in the limit ^ — > 0. Eq.( |3.28| ) is nothing but the sum of 
two quantities 

= V{m;0 + V{m;^,oo,fi) , (C.l) 

where V{m,^) is given by ( p.27| ), and V{m; ^, oo, fi) is given by RHS of ( |3.22| ). According 
to the result of Appendix B, we have 

tr[l]2v/^ fi^ 



S^^"^'^'-'^)=(4.)-/^r(i^)D(D-i)- 



(C.2) 



Let us consider the limit of another quantity V{m; ^). In the situation of ^ —^0, three 



terms of zeta function on RHS of (|3.27|) only survive 
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and we adopt the following continuum limit as ^ — ^ 0, 
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The limit of ( |C.3| ) is therefore led to the following 



lim V(m: 
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Combining ( |C.2| ) and (|C.5| ), the limit of turns out to be 
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This is equivalent to ( p.31| ), and the proof ends. 
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D The explicit forms of hi^o-) and 



We write down explicit expressions for 12(0") and h{cy)i which are defined by 
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And similarly, 13(a) is composed of: 
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Our notation for the Lerch transcendental ^{z\ s, a) is 
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and a Lerch transformation is then expressed by 
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Taking a = 1, we have 
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